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RÉSUMÉ 


L'article comprend deux parties. 

Dans une première partie, les auteurs analysent les bases physiques et mathé- 
matiques des méthodes de mesure de la conductivité capillaire à partir de la 
vitesse de drainage d’un échantillon de sol soumis à une brusque augmentation 
de tension. Ils comparent les résultats expérimentaux aux conséquences théo- 
riques des lois utilisées. Ils aboutissent à des discordances importantes justifiant 
le rejet des hypothèses. 

Dans une seconde partie, les auteurs émettent l'hypothèse que le potentiel de 
l’eau liée au sol, au cours du phénomène envisagé, varie non seulement en fonction 
de la teneur en eau mais également en fonction des liaisons dynamiques, variables 
au cours du temps, qui s'appliquent à l’ensemble du système. 

Ces nouvelles hypothèses sont en meilleure concordance avec les résultats 
expérimentaux. 

Elles ne sont pas en opposition avec la loi de Darcy. 
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NOTATIONS UTILISÉES 


volume de l'échantillon en em. 


: section en cm. 

: hauteur en cm. 

: temps. 

: pression. 

: variation de pression. 

: teneur en humidité du sol exprimée en volume d'eau par volume 


apparent de sol. 


: volume d'air par unité de volume apparent de sol. 

: débitvolumétrique dans la directions, par unité de section, normale à s. 
: coefficient de conductivité hydraulique ou capillaire. 

: altitude de la section d'écoulement. 

: direction et axe de l'écoulement. 

: potentiel de capillarité, exprimé en hauteur d'eau, des points d’alti- 


tude z. 


: coefficient de diffusivité capillaire. 


base des logarithmes népériens. 


: chute de potentiel piézométrique au travers de la membrane ou de la 


plaque. 


: potentiel piézométrique = somme de l'altitude et de la hauteur 


piézométrique. 


: épaisseur de la membrane ou plaque poreuse. 
: conductivité hydraulique (coefficient de Darcy) de la membrane ou 


de la plaque poreuse saturée. 


: impédance de la membrane. 


: tension superficielle. 

: angle de contact. 

: quantité d'eau écoulée au temps t. 

: quantité totale d'eau écoulée. 

: surface normale aux lignes de courant, par unité de volume au point 


d'altitude z. 


: volume de la partie solide, rapportée à l'unité de volume apparent 


de sol. 


: rayon de courbure principal des ménisques. 
: rayon des cavités gazeuses correspondant à la formule (45). 
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Tre PARTIE 


LES BASES DE LA MÉTHODE GARDNER 


La mesure de la conductivité hydraulique ou capillaire des sols 
s'obtient, en milieux non saturés, à partir de plusieurs méthodes. 
Parmi celles-ci, une des méthodes les plus courantes, actuellement, 
est celle imaginée par GARDNER (3) et adaptée par BuTIJN et WessE- 
LING (1) pour les teneurs en humidité élevées. Ces auteurs considérent 
un échantillon cylindrique de sol, de volume V, de section A et de 
hauteur L. Cet échantillon repose sur une membrane (ou une plaque) 
poreuse, horizontale, perméable à Peau et imperméable à Pair. 
Gardner utilise la membrane de l'appareil à pression de RicrarDs (7). 
Il provoque à l'instant t = 0 une augmentation AP de la pression 
initiale P; règnant à l’intérieur de la chambre. Avant cet instant, l’échan- 
tillon humide est en équilibre à la pression P, et sa teneur en humidité 


est uniforme et égale à 6. 


Schéma de l'appareil de Butijn-Wesseling 


anneau cylindrique 
de section À 


plaque 
poreuse semi- 
perméable 


BA) 


tube calibré 


goutte de mercure 


D CRAF(IRS..A ) 
Gembloux 
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L'augmentation de pression de lair, à l’intérieur de la chambre, 
provoque un écoulement d'une partie de l’eau de l'échantillon, au 
travers de la membrane poreuse. L’écoulement s'arrête lorsque la 
teneur en humidité de l’échantillon atteint une nouvelle teneur en 
humidité 0, en équilibre avec la pression P; + AP. 

Bruns et WESsELING utilisent une plaque poreuse, analogue à celle 
d’un tensiométre (fig. 1). La teneur en humidité initiale de l'échantillon 
de sol est homogène, égale à 6,, en équilibre avec la dépression initiale D. 

On provoque à l'instant t = 0, une augmentation AP de la tension 
initiale P. L'augmentation de tension entraîne un écoulement qui cesse 
lorsque la teneur en humidité de l’échantillon atteint, partout, une 
nouvelle teneur en humidité 0,, en équilibre avec la tension P; + AP = 
P,. Cet écoulement est observé à partir du déplacement d’une goutte de 
mercure dans le tube calibré. Les lois qui régissent les écoulements dans 
les 2 types d’expériences (GARDNER et BUTIJN) sont, d’après ces auteurs, 
analogues au point de vue formel. Elles permettent le calcul du coeffi- 
cient de conductivité capillaire à partir des raisonnements ci-dessous. 


MISE EN ÉQUATION 


On admet que le mouvement unidimensionnel de l’eau, dans un 
milieu filtrant non saturé est régi par la loi généralisée de Darcy : 


q=—K°_@—y) (1) 
Os 
dans laquelle : 
q — débit volumétrique dans la direction s, par unité de section, 
normale a s. 
K= coefficient de conductivité hydraulique ou capillaire. 
z = altitude de la section d'écoulement. 
s = direction et axe de l’écoulement. 
y = potentiel de capillarité des points d’altitude z exprimé en 


hauteur d’eau. 


On admet également l'équation de continuité d'un écoulement uni- 
dimensionnel : 


60 Ô 
a" ie =% 2) 
dans laquelle : 
t = temps 
0 = teneur en humidité du sol, exprimée en volume d'eau par 


volume apparent de sol. 


396 Bull. Inst. agron. et Stat. Rech. Gembloux TOME XXXI, N° 3 


Dans un intervalle d'humidité 0, — 0; suffisamment faible, on admet 
la constance du coefficient K. 

Dans ces conditions, la combinaison des équations (1) et (2) aboutit à 
la relation : 


LK (3) 


GARDNER, ainsi que Burg et WESSELING, admettent que, durant 
tout l'écoulement, le potentiel y n’est fonction que de la teneur en 
humidité 0 et que la relation reliant w et 0 est connue. Cette relation 
est la relation classique définissant le pF, puisqu'on a (8) : 


pF = logi 


lorsque y est exprimé en cm. 


Ils posent : 
dy 
D=—K EK (4) 


D étant appelé le coefficient de diffusivité capillaire. 
L’équation (4) se transforme dès lors en l'équation de la chaleur (ou 

de la diffusion) 
òG D 0 


ot os? 


(5) 


équation aux dérivées partielles qu’il faut résoudre en fonction des 
conditions initiales et des conditions aux limites. 


MEMBRANE (OU PLAQUE) POREUSE DE PERMEABILITE INFINIE 


Lorsque la perméabilité de la membrane (ou de la plaque) est négli- 
geable, les conditions initiales et aux limites des expériences ci-dessus 
peuvent s'écrire : 

Dis, o) = 6, pour o < s < L 


Do, t) = 0, pour t > o (6) 
60 0 
ds W 


s=L 


en prenant pour axe des s, l'axe vertical orienté vers le haut et ayant 


son origine dans le plan séparant la membrane de l’échantillon (axe 
des z). 
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La solution de l’équation (5), satisfaisant aux conditions (6) est classi- 
que (cf. CARSLAW et JAEGER) (2, p. 79). 


0-0, = 4 nm “FE 7 
0,-0 lag DR sin JL c (7) 


=y 


f n= 


On peut en déduire la loi du flux de drainage de Péchantillon à partir 
des relations (1) et (4). 
00 


òs R (8) 


qs-0 = — 
On admet généralement que la valeur de K est négligeable devant 


id 
le terme D ss et l’on calcule 
o 


òn 


Os (8a) 


Han — 


Dans ces conditions, la quantité d'eau Q,, qui quitte l'échantillon 


entre l'instant t = o et l'instant t, est fournie par : 
wi 
RI 
= | T (E dt (9) 
Qı Jo (3s a 
c'est-à-dire : 

g DN E 

Q e 4L \ (10) 


dans laquelle Qt représente la quantité totale d’eau écoulée et est égal à 
Q; = AL. (0; —0;) 
La série (10) converge rapidement et pour t suffisamment grand, on 
peut se contenter d'utiliser le premier terme de la série et écrire : 
0 _ 8 + (11) 
J =o ==; e 
Q 7 
équation qui permet de calculer facilement D, à partir d’un alignement 


j e t A x 
des points représentant | — Q en fonction de t dans un diagramme 
f 


semi-logarithmique. 
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Pour des valeurs de t proches de 0 on peut recourir à une autre 
solution approchée du problème. Elle consiste à admettre, qu’en début 
d'écoulement, celui-ci correspond à l'écoulement dans un milieu semi- 
infini satisfaisant à l'équation fondamentale (5) et aux conditions 
initiales et aux limites 


O(s, 0) = 0, pour 0 < s (6a) 
0(O, t) = 4, pour 0 < t 


La solution (2, p. 41) est donnée par : 


9 — 0 
A = èïf u (12) 
moyennant : 

> 13 

u — 
2VDt (23) 

2 CR a 
erf u = = e" du 
VS Jo 


La quantité d’eau écoulée, jusqu’à l'instant t, calculée à partir des 
relations 9 et 12 aboutit à: 


2 = 
Q= RE He V Di (14) 


Pour comparer la loi (10) aux approximations 11 et 14, nous poserons : 
Dt 
Et (15) 


Les lois (10) (11) et (14) s’écrivent respectivement : 


Qt 8 = Le 
1 D = a =, Gg t =1— F(x) (10a) 
wx 
e = 
SET ae? = 1—F,(x) (la) 
Q: 2 = 
SE = F,(x) (14a) 


Le tableau I indique les valeurs respectives de F, (x), F, (x) et F, (x). 
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TABLEAU I 


x I-F, (x) I-F, (x) I-F; (x) 
0 1,0000 1,0000 
0,05 0.7477 0,7477 
0.10 0,6432 0,6432 
0.15 0,5631 0,5630 
0.20 0,4960 0,4954 
0.3 0.3868 0,3867 0,3820 
0,4 0.3021 0.3021 0,2863 
0.5 0,2361 0,2361 0,2021 
0.6 0,1844 0,1844 0.1260 
0,7 0,1441 0.1441 0,0559 
0.8 0,1126 0.1126 —0,0092 
0.9 0,08798 0.08798 — 
1.0 0,06874 0.06874 — 
œ 0 0 — œ 


Ce tableau démontre qu'on peut toujours recourir à une des formules 
:implifiées (11) ou (14) pour autant que l'on respecte les conditions 
suivantes : 


Dt er 
si Le est inférieur ou égal à 0,2: formule 14 
. Dt BS n 
si Co est supérieur à 0,2 : formule 11 


CALCUL DE LA DIFFUSIVITÉ CAPILLAIRE 


Les formules 11 et 14 permettent, quel que soit t, un calcul aisé du 
coefficient de conductivité capillaire qui s'obtient à partir de ces rela- 
tions que lon écrira sous la forme : 


pour a < 0,2 (faibles valeurs de t) 
= ys (QU L. 
Be eles (14b) 
Pour 2 > 0,2 
Qt 
DR 412 A log, (1 QQ (11b) 
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= D 
MEMBRANE (OU PLAQUE) POREUSE DE PERMBABILITÉ FINIE 


Cet aspect du problème a été envisagé par MILLER et ELRICK (4) et 
par RijremA (5), à partir des mêmes hypothèses de calcul. 

Ces auteurs font remarquer qu’il y a lieu de modifier une des condi- 
tions aux limites contenues dans les relations (6) afin de tenir compte 
de la perméabilité finie de la membrane ou de la plaque. 

La deuxième condition des relations (6) sera modifiée. 

On exprime l'égalité des flux d'écoulement au travers de la section de 
contact échantillon-membrane à partir de la combinaison de la rela- 
tion (8a) et de la relation exprimant l'application de la loi de Darcy, 
(en milieu saturé) à la membrane. 


On écrira donc : 
- =- — Aso (16) 


dans laquelle : 


lo = chute de potentiel piézométrique au travers de la membrane 


(plaque) 

ø = potentiel piézométrique, somme de l'altitude et de la hauteur 
piézométrique. 

Lm = épaisseur de la membrane poreuse (plaque) 

K,, = conductivité hydraulique (coefficient de Darcy) de la mem- 


brane (ou de la plaque poreuse) saturée. 


Dès lors, les nouvelles conditions initiales et aux limites s’écrivent : 


0.) = 0; pour o <s < L 

ag Ky 

Bs. = LE ` 3 (9,4 — 05) (17) 
50 

ei = Ü 

Os, 


s= L 


K 


m 


Le rapport a été appelé « impédance» de la membrane (4) (5). 


Il faut remarquer qu'en écrivant la deuxième condition de la relation 
(17) MILLER et Evrick, ainsi que Riyrema, supposent implicitement 
que la relation (16) puisse s'écrire indifféremment : 


oy Ys=o — Wr 00 K Ps=0 — Pr 


ES = — 0; _ 0 
ds d as fr dë L. Vg ht (Bu — 0) 


s 
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c’est-à-dire : 


K Z Se (9 9 aa 05) (16a) 


qui est la deuxième condition de la relation 17. 

Nous reviendrons ultérieurement sur la validité de ces hypothèses. 
A partir des conditions (17) l'équation fondamentale (5) a pour solution 
(cf. (4) (5) et aussi (2 p. 101). 


_ an Dt 
Sy], Le 
o A , e cos [ (1 all (18) 
0; 0; ` g= 1,2. ay sin Wei (a kb: cosec?a,) 
_ mé Di 
LE 
RC 8 D 
Q: n=1,2.. A (a + cosec?a,) 
moyennant les notations : 
K La 
a= Ç K, (20) 
a, = n®™e racine de Péquation 
aa, tga, = 1 (21) 
Si Pon pose, en outre : 
vDt 
RATES" (22 
aL % ) 
l'équation (19) s'écrit : 
Qt = — du ay? 
1 = = y A e (19a) 
Qr a nf 


les coefficients a, étant définis par la relation (21) et les coefficients A, 
étant égaux a 

2 2 
Tis TL: (23) 


"` a (a + cosec? a) aa? (a2 + la + 1/a?) 


A 


et étant fonction de n et de a. 


Les valeurs des coefficients A, correspondant à n variant de l à 6 
peuvent être calculées à partir des solutions a, données par CARSLAW 
et JAEGER et de la relation 23. 
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Comme précédemment (formule 10), on constate que la série (19a) est 
rapidement convergente pour autant que t soit suffisant. 

Si l’on se fixe une précision de 1%, on peut, dans l’équation (19a), 
ne retenir que le premier terme, pour autant que 


x = ap = 202 (24) 


Le tableau (II) le démontre à partir du calcul numérique des valeurs 


t e . 
de 1 — Qu en fonction de a et de y et des relations : 


Qr 
eo 2 
Lass É =F | A, c ™* = 1—F, (a,x) (19b) 
f n=1,2... 
ï — E =A,e %* = 1—F, (a,x) (25) 
f 
a fonction F, correspondant à: 
=A 
Qu S e a 


1 


Q: 


a, étant la première solution de équation (21). 
d D q 


les Fi e= (25a) 


a (a + cosec?a) 


Pour les faibles valeurs de t, on peut écrire les conditions aux limites 
et initiales (17) sous la forme 


(so) = 0, Douro < s 
OG 1 

=> 5 | ech 7 
ER T Cot) (17a) 


La solution mathématique est indiquée par CARSLAW et JAEGER (2, 
p. 51). 
= S _ 4 gal tan * vDi 


erf — 
0; — 6; 2 VDt 


EE 


moyennant la notation (13) et la notation complémentaire (2) 
^ 


2 e 
-= | ST du 
u 


erfc u = 1 — erf u 


` Vx. 
La quantité d'eau écoulée jusqu'à l'instant t calculée à partir des 
relations (9) et (26) aboutit a: 
Dt — 
TS Ji 2 = 
Z [et ete VD 4 E vp — 1] = 27) 
a 7 
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Le tableau II montre que cette formule est valable, avec une précision 
supérieure à 1 %, pour autant que 


La figure 2 représente graphiquement la variation de 1— sen fonction 
$ 


de x, calculée à partir de la relation générale (19b). 
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TABLEAU Il 


a LP, 1-F, I-F, 
x =:0,1 
0 0,1432 0,6334 0,6432 
0,05 | 0,6888 0.6767 0,6887 
0,10 0,7261 0,7128 0,7262 
0,20 0,7814 0,7683 0,7814 
0,50 0,8663 0,8580 0,8664 
1,00 0,9197 0,9158 0,9196 
2,00 0,955: 0,9539 0,9554 
5,00 0,9809 0,9807 0,9809 
10,00 0,9902 0,9901 0,9902 
x = 0,2 
0 0,4960 0,4949 0,4954 
0,05 0,5426 0,5410 0,5422 
0,10 0,5832 0,5812 0,5830 
0,20 0,6491 0,6465 0,6489 
0,50: 0,7663 0,7640 0,7662 
1,00 0,8517 0,8505 0,8516 
2,00 0,91485 0,9141 0,9146 
5,00 0,9625 0,9624 0,9625 
10,00 0,9806 0,9806 0,9806 
x= Q05 
0 0,3868 0,3867 0,3820 
0,05 0,4327 0,4325 0,4294 
0.10 0,4742 0,4739 0,4718 
0,20 0,5445 0,5440 0,5431 
0,50 0,6810 0,6804 0.6805 
1,00 0,7902 0,7898 0,7899 
2,00 0,8760 0,8759 0,8761 
5,00 0,9446 0,9446 0,9446 
10,00 0,9713 0,9712 0,9712 


Cette relation s'écrit, en utilisant la notation 


Qt 
aQ; 


2 
= [ d erfc y + 


2y 


VT 


(22) 


dé 


(27a) 


Nous avons posé, par analogie aux relations 14a, 19b et 25, 


= S =1=a [er 


2 
erfc (xl? a-1) + 


232 =1] KH 


aya 
1—F, (ax) (28) 
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DETERMINATION DE a, 


Diverses méthodes ont été proposées. 

La méthode la plus séduisante, sur le plan mathématique, nous 
paraît être celle proposée par P. V. Riyrema. 

Il remarque que, dans un diagramme semi-logarithmique, les droites 
qui représentent les équations (11) et (25) coupent l'axe t = 0 en un 

2 
aï (a + coscc? a), 
c’est-à-dire en un point qui est indépendant de D et ne dépend que du 
facteur a. Il n’est dès lors pas difficile théoriquement, à partir d’un 
graphique expérimental et d’un tableau de correspondance A, (a) et A, 
(a) (cf. tableau IIJ), de calculer œ et a, sans passer par la mesure de 
l'impédance de la membrane. 

En pratique, le faible intervalle de variation de A, (cf. tableau III, 
A, variant de 0,81 à 1,00 pour « a» variant de 0 à P œ) et l'allure des 
courbes expérimentales (cf. figure 3) rendent utopique l'utilisation de 
cette méthode. 

MILLER et ELRICK (4) remarquent que l'on peut écrire, pour une 


valeur suffisamment élevée de t (c'est-à-dire en combinant les relations 
21 ét 25a). 


point d’ordonnée égale respectivement à ou à 


a 


a 


tga, 


KAL D LL,, 4-9; 99 
~ 84 TIK, ` L? Ky sn e 


m 


s do A e , 9; 0 
si l'on admet que la grandeur => (relation 4) est assimBable à == 
du Vi Yr 


Dès lors, ils remarquent que la valeur : 
D LL 9-9; 
L? Ka Pv 


m 


peut être calculée à partir d'un diagramme à échelle semi-logarithmique 


sur lequel on a reporté log, (1 — en fonction de t. Le coefficient 


Q: 

Of) 
angulaire de la droite exprimant l’équation (25a) fournit la valeur de 
ai D 
L: 
mesurer. 


; les autres termes du second membre de l'équation (29) peuvent se 


a = 
Ayant calculé la valeur = , Muer et Evrick calculent la valeur 
a 


de D à partir d'un procédé graphique. On peut également procéder par 
voie analytique ainsi que nous le démontrons ci-dessous. 
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Connaissant, par la méthode de Miter et ELRICK, la valeur de = 
1 


on déduit a, par le tableau de correspondance que nous avons calculé 
(tableau III). Connaissant œ, et la valeur précédemment calculée de 
d 


, on déduit la valeur de D. 
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Tasreau III 


a 2 


A ee! 
tga, = a2(a+ cosec*a,) 


tga, | aiga, 


0 | 0 1,000 1,0000 
0,0100 0,0001 8 1,000 1,0000 
0,0200 0,0004 R 1,000 1,0000 
0,0300 0,0009 À 1,000 1,0000 
0,0400 0,0016 8 1,000 1,0000 
0,0601 0,0036 0,999 1,0000 
0,0802 0.0064 á 0,998 1,0000 
0,1003 0,0100 0,997 1,0006 
0.1511 0,0227 0,992 1,0600 
0,2027 0,0405 2 5 0,987 1,0000 
0,2553 0.0638 5.66 0,979 0.9999 
0,3093 0,0928 i 0,970 0,9998 
0,4228 0,1691 5; 0,946 0,9994 
0,5463 0.2731 .662 0.915 0,9985 
0,6841 0.4105 2,436 0.877 0,9969 
0,8423 0.5896 .696 0.831 0,9942 
1.0296 0,8237 # 0,777 0.9898 
1.2602 1.1341 H 0,714 0.9831 
1,5574 1,5574 6 0,642 0.9735 
1.7433 1,8305 j 0,602 0,9673 
1,9648 2.1612 0,560 0,9600 
2,2345 2,5699 0,515 0,9515 
2,5721 3,0866 0,466 0,9415 
3,0096 3.7620 267 0,415 0,9300 
3.6021 4,6827 ` 0,361 0.9169 
4.4552 6.0145 ,166 0,303 0,9020 
5,79798 8.1170 0,1232 0,241 0,8850 
8,23819 11,9452 0,08372 0,176 0,8680 
8,98861 13,1234 0,07620 0,162 0,8619 
9,88747 14,534 0,06880 0,149 0,8578 
10,98348 16,255 0,06152 0.135 0,8535 
12,34996 18.401 0,05434 0,121 0,8492 
14,101 21,152 0,04728 0.106 0,8447 
16,428 24,806 0,04031 0,092 0,8402 
19.670 29,898 0,03345 0,077 0,8355 
24,498 37,483 0,02668 0,062 0.8308 
32,461 49,990 0,02000 0,047 0.8260 
48,078 74,522 0,01342 0,032 0,8211 
92,620 114,488 0,00873 0,017 0,8146 
1255,766 1971,553 0,000507 0,001 0,8110 


co Ka 0 0 0,8106 
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RÉSULTATS EXPÉRIMENTAUX 


Nous avons voulu vérifier l’utilisation des théories ci-dessus, afin de 
déterminer expérimentalement le coefficient de conductivité capillaire 
de différents sols. 

En raison des difficultés d'interprétation des résultats expérimentaux 
obtenus, nous avons finalement étudié un milieu aussi simple que pos- 
avons eu recours à un processus expérimental également 


sible et nous 
aussi dépouillé que possible, 

Nous avons utilisé un milieu composé de billes de verre, de diamètres 
assez uniformes, compris entre 125 et 250 y. 

Nous avons effectué les mesures à partir du schéma imaginé par 
Burin et WeEssELING (1) en remplaçant la goutte de mercure par une 
bulle d’air dont le déplacement est contrôlé en fonction du temps 


4). 


(cf. fig 


PRET anneau cylindrique 
Vue de l'appareil de section A 


j/aque poreuse 
perméable 


tube calibré 
` Fig. 4 


Laf courbe de y en fonction de 0 de ce milieu, pour des tensions va- 
riant de O à 160 cm (c'est-à-dire entre pF 0 et pF 2,2) est ren- 
seignée figure 5. 
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On remarquera que les mesures expérimentales se rapportent unique- 
ment à des teneurs en humidité élevées (pF < 2,2). Le diagramme de 
la figure 5 résulte en moyenne de 7'mesures par point, la dispersion des 
résultats étant indiquée sur le diagramme par la représentation de 
l'intervalle de variation des mesures. 


H 


L'apparcillage utilisé est le même que celui que nous employons pour 
la mesure de la conductivité capillaire. 

Pour compléter ces mesures de la courbe y (0), obtenues à partir 
d'échantillons de 1 em. de hauteur, nous avons, en outre, déterminé les 
profils hydriques dans des échantillons de 12 cm. de hauteur soumis à 
diverses tensions y à la base de l'échantillon. 

Les profils hydriques ainsi obtenus sont représentés à la fig. 6. 

On peut constater que, dans la zone de forte variation de y en fonction 
de 6 (c’est-à-dire pour y compris entre 40 et 80 cm), il apparaît un 
profil hydrique différencié et on peut y admettre, à l'équilibre, la relation 
hydrostatique des pressions 


(cf. fig. 6). 
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il 


$ Z 6 7 8 9 10 n Iz VE 16 17 18 19 20 À 2 23 24 26 
Teneur en humidité en ‘lo du poids sec 


Les résultats des mesures d'écoulement ayant servi à la détermina- 
tion du coefficient de diffusivité capillaire sont illustrés par les quelques 


courbes de la figure 3. L'ensemble des mesures est résumé au tableau 
IV. 


On y trouve successivement, classés par colonne : 


— la numérotation de l'essai. 
— la hauteur de l'échantillon L en cm. 
— les valeurs des potentiels initial et final, 
y et y, en cm 
— le poids de la matière sèche P, en gr. 
— le poids spécifique apparent de la matière sèche en gr em? 
— la quantité d’eau écoulée Q, en cm? 


— le coefficient angulaire de la droite représentant loge (1 -3 
ï 

aD 

L2 


en fonction de t c’est-à-dire : 
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a à a 
— le coefficient ` 
tga, 


q 


— le rapport Sa entre la vitesse initiale d'écoulement et la 
max 


vitesse maximum théorique qmax calculée par la relation (31) 
— le coefficient D calculé par la relation (11 b) 
— le coefficient D U calculé par la relation (43) 
— le coefficient D calculé par les relations (43) et (45). 


On remarquera que les mesures ont porté sur deux types d’échan- 
tillons de hauteurs respectives égales à 1 et 3 cm. 

La section de tous les échantillons était constante. On a utilisé des 
tubes calibrés en acier de section égale à 17,94 cm. 

Nous n’avons pas reporté les mesures des écoulements inférieurs à 
0,75 % du poids sec. Ces écoulements correspondaient le plus souvent 
à des potentiels inférieurs à 50 cm ou supérieurs à 100 cm. La différence 
vc est toujours égale à 10 cm. pour tous les écoulements observés. 


Ile PARTIE 


CRITIQUE DES HYPOTHÈSES DE CALCUL ET INTERPRÉTATION DES 
RÉSULTATS EXPÉRIMENTAUX 


Considérons le cas d’une plaque de perméabilité finie et envisageons 

successivement : 

— la simplification apportée par la formule (8a) par rapport à la 
formule (8), c’est-à-dire l'importance comparative du phénomène 
de percolation par rapport au phénomène de diffusion. 

— les conséquences qu’impliquent la 2° condition des relations (17) 

— les possibilités d'interprétation de la loi (1). 


PERCOLATION ET DIFFUSION 


Si on applique la loi (1) en toute rigueur (formule 8), on admet que 
l'écoulement cesse dans l'échantillon lorsque l'influence des forces de 
gravité est équilibrée par l'influence inverse d’un gradient des forces de 
capillarité, 

Après arrêt de l'expérience, un profil hydrique plus ou moins diffé- 
rencié indique donc l'importance relative des phénomènes de perco- 
lation par rapport aux phénomènes de diffusion. 


TABLEAU TV 


D T E =Ë = 
li 1 lh] 4 dr a 
ERR Sal wël $ A ER 5 z 
S | ES | 2 3 v| & ï k M 53 F 3 CS 
v z Ki Boe z 5 A a= > Se 
8 [ëg] E. |Ë E| ï| aa 3 iS # BST 6 =* 
2 ERT a E z B ë É F 8 £ 5 E i =| Ë É a s Ê 3 
Ë  385|/ 5/ 285 38. %⁄| 38 | a e "| sË z 5 Ë 
5 2/5] & se | 55 ERS SI slg al E 2 SE 
Z | 23 | És és fe) a8 | Gë vi ZS as A as 
1 2 3 4 5 6 7 8 9 10 11 12 13 
1 1 50 60 | 29,68 | 1,65 1,180 83,1 1,75 0,400 33,7 32,6 0,43 
60 70 » 1,970 113,8 4,05 0,555 46,1 37,0 0,39 
70 80 » 1,110 255,7 5,22 0,670 103,6 92,5 0,85 
80 90 » 0,495 519,3 4,74 0,593 210,5 185,0 1,56 
90 100 » 0,309 692,4 3.96 0,357 280,6 277,0 2,17 
2 1 50 60 | 30,05 | 1,70 | 1460 61,6 1,67 0,705 25,0 31,0 0,40 
60 70 » 2400 78,5 3.48 0,765 31,8 24,2 0,25 
70 80 » | 1,180 221,6 4,95 0,892 89,8 79,2 0,72 
80 90 » 0,505 401,6 3,70 0,595 162,8 158,0 1,32 
90 | 100 » 0,337 426,9 2,65 0.537 173,0 201,5 1,58 
3 1 50 60 | 2845 | 1,59 | 1,430 63,0 1,67 0,473 25,5 27,6 0,36 
60 70 » 1,810 117,7 3.90 0,640 47,7 41,2 0,43 
70 80 » 1,130 267,7 5,52 0,596 108,5 79,2 0,72 
80 90 » 0.475 396,9 3,48 0,430 160,9 158,0 1,32 
90 | 100 » 0,235 461,6 201 0,179 187,1 185,0 1,46 
4 1 50 60 | 29,09 | 1,62 1,240 448 1,02 0,330 18,2 19,2 0,26 
60 70 » 2,030 85,6 3,24 0,680 34,7 33,6 0,36 
70 80 » 1,420 188,8 5,00 0,670 765 65,1 0,59 
80 90 » 
90 100 » 0,500 238,2 2,32 0,417 96,5 158,0 1,25 
5 1 50 60 | 26,10 | 1,45 | 0,900 13,9 0,38 0,165 5,6 17,3 0,23 
60 70 » 1,980 25,2 1,38 0,283 10.2 18,5 0,20 
70 80 » 1,170 186,7 6,05 0,740 75,7 50,4 0,46 
80 90 » 0,252 280,2 1.96 0,110 113,6 79,2 0,61 
90 100 » 0,057 179,6 0,28 


bag 


XNO]QUIE) “YHA “WIS 12 NIT "pul "ng 


$ oN ‘IXXX ANOL 


1 2 3 4 5 6 7 8 9 10 11 12 13 
6 1 50 60 | 27,74 | 1,55 0,423 116,1 0,91 0,470 47,1 85,5 1,14 
60 70 » 2,320 70.6 3,03 0,680 30,8 26,3 0,29 

70 80 » | 1,670 190.9 5.90 0,745 774 65,0 0,59 

80 90 » | 0,472 405,3 3,53 0,307 164,3 158,0 1,32 

90 100 » 048,3 503,4 3,23 0,356 204,0 222,0 174 

7 1 50 60 | 27,04 | 1,52 0146 722,4 1,96 0,375 27,2 = = 
60 70 » | 3320 67.2 4.12 0,854 102,0 23,6 0,26 

70 80 » 1095 251,8 5,10 0.600 177,3 85,0 0,78 

80 90 » 0.510 437,4 4,11 0,529 292,8 185,0 1,56 

90 100 » 0.258 722,4 3.45 0,387 227,2 370,0 2,88 

8 3 50 60 | 85,17 | 1,58 1,810 62,3 2,10 0,845 227.2 101,0 1,34 
60 70 » 4,850 44,3 3,96 0,768 161,6 46,2 0,49 

70 80 » 4.000 68,1 5,03 0,820 248,4 69,3 0,63 

80 90 0,124 662,4 1,53 0,132 = = = 

90 100 » 0,062 920,0 1,03 0,157 f = = 

9 3 50 60 | 85,37 | 1,59 | 3,040 17,2 1,44 1,030 62,9 20,7 0,27 
60 70 » 4350 24,2 2,80 1,161 88.3 34,0 0,31 

70 80 » | 5400 41,5 6,18 0,945 1514 41,8 0,39 

80 90 » | 1,770 97,5 4,58 0,515 355,6 83,0 0,71 

90 100 » 0,825 106,0 241 0,335 386,6 128,0 1,01 

10 3 50 60 | 84,44 | 1,58 | 4,260 19,7 0,375 719 22,2 0,30 
60 70 » 3,330 31,0 0,683 113,1 37,0 0,39 

70 80 » 2,660 75,0 0410 273,6 75,5 0,69 

| 80 90 » 1,520 108,0 0,364 393,9 119,0 0,99 

90 100 » 0,790 129,5 0,297 4724 166,0 1,31 

11 3 50 60 | 82,63 | 1,54 2,680 344 2,55 0,750 125,5 43,5 0,58 
60 70 » 5,550 29,3 4,50 0,725 106,9 37,0 0,39 

70 80 » 3400 66,5 6,23 0,374 242,6 66,5 0,60 

80 90 » 1,530 119,5 5,00 0,465 435,9 119,0 0,99 

90 100 » 0,900 117,0 2,90 0,253 426.8 119,0 0,93 

12 3 50 60 | 87,19 | 1,51 4,350 29,8 3,57 0,845 108,7 39,5 0,53 
60 70 » 4,550 43,8 5,50 0,810 159,8 52,0 0,55 
70 80 » 3,440 67,2 6,38 0,523 245,1 72,0 0,66 
80 90 » | 1,085 74,8 2,24 0,397 272,8 92,0 0,77 

90 100 » 0,186 449,1 2,30 0,238 = — 


SE —_—]Â © Ù Ù — 


Remarque : 


q exprime le rapport entre les vitesses d'écoulement au début de l'expérience (au cours de la première minute) 


max et la vitesse maximum d'écoulement calculée par la relation (31). 


sjos sap auvpdva gnoympuoy = "V ‘ZINI 79 "T “INIS 


ER 
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Nous avons vu que, dans la zone des écoulements importants, cette 
influence peut être sensible (cf. fig. 6). 

On ne peut la négliger, comme dans l’utilisation des formules établies 
ci-dessus, que si l'échantillon n’est pas de hauteur L trop importante 
par rapport à la tension v. 

Dans la gamme des résultats expérimentaux envisagés (L variant de 
1 à 3 cm), l'influence de forces de gravité paraît négligeable. 


POTENTIEL DE L'EAU DANS LA PLAQUE 


On peut se demander si la deuxième condition de l'équation (17) telle 
qu'elle est exprimée par l'expression (16a), correspond à la réalité. Cette 
équation résulte de l'application de l'équation de continuité des débits 
au travers de la section de contact entre l'échantillon et la plaque 
poreuse. Elle est basée sur deux hypothèses différentes concernant 
respectivement l’état énergétique de l’eau dans la plaque et dans l’échan- 
tillon, à la surface de contact. 

Elle suppose, en premier lieu, que toutes les particules d’eau à la face 
supérieure de la plaque ont un potentiel ø, calculable par la relation 
suivante : 


a=— +| g” (30) 
m 
qı représentant le débit volumétrique traversant la plaque, par unité 
de surface à l'instant t. 
Elle suppose, en outre, qu'en tous ces points (face supérieure de la 
plaque) le potentiel ø, varie de —y, à _ au cours de l'écoulement. 
En conséquence, la vitesse maximum d'écoulement au travers de la 
plaque est donnée par la relation 


K 
gen = (Vr V) T (30a) 
m 


On peut, cependant, tenir un raisonnement différent et aboutir à 
une formule différente pour Oe, formule en contradiction avec la 
relation (30a). 

Avant le début de l'expérience, la plaque supporte une tension corres- 
pondant à y, Elle est saturée et ne peut être le siège d'un gradient de 
pression, puisqu'il n'y a pas de mouvement. Dès lors, les ménisques 
supérieurs, à la face de contact avec l'échantillon doivent avoir un 
rayon r; correspondant à y, c'est-à-dire être en équilibre avec y. 
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Dès lors, si nous admettons qu’au niveau de la plaque on rencontre 
des ménisques de rayon r; tels que: 


26 cosb 
q= 
Vi 
dans laquelle : 


a = tension superficielle 
b — angle de contact 


la mise sous tension y, va provoquer l'apparition instantanée de nou- 
veaux ménisques r, en équilibre avec w, et, lorsque l'écoulement se 
produira, on peut émettre l'hypothèse que l'apport d’eau provoque la 
disparition momentanée de ces ménisques et que le gradient de pression, 
au travers de la plaque, est initialement ze Ly, 

Dans ces conditions, la vitesse maximum devient, pour les expériences 
étudiées : 


=7= M (31) 
max SS A ) 
formule qui est en contradiction avec la relation (30a). 


Or nos résultats expérimentaux sont en contradiction formelle avec 
la loi (30a) base de la méthode de calcul proposée par MILLER et ELRICK, 
ainsi que, par Rijrema (cf. tableau IV). 

En d'autres termes, la deuxième condition aux limites des équa- 
tions (17) ne correspond pas à la réalité. Il cn est de même, en consé- 
quence, des formules d'intégration utilisées par MILLER et ELRICK ou 
par RIJTEMA. 

On trouvera une preuve complémentaire de notre assertion dans le 
fait que, dans nos résultats d'observation, le rapport a (gu calculé à 
partir de la relation (29) est souvent supérieur à 1 (cf. tableau IV) 
démontrant, par là même, l'impossibilité de déterminer œ, la valeur 


d a 2 à e 
théorique de = devant toujours être comprise entre 0 et 1 (cf. ta- 
a: 


1 
beau III). 
En conclusion, l'influence de la perméabilité de la membrane sur le 
mode d’écoulement est nettement inférieure à Vinfluence calculée, 
explicitée par les équations (19) ou (26). La différence est d'autant plus 


importante que le rapport © — 1 est petit. 
t 


Dans ces conditions, si les hypothéses de calcul (1) et (2) sont exactes, 
il est préférable de recourir systématiquement, aux équations (11) ou 


(14). 
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VALIDITÉ DES HYPOTHESES EXPLICITÉES PAR LES LOIS (1) ET (2) 


Dans l'interprétation proposée par GARDNER et par BUTIJN et WESSE- 
LING, on admet que la vitesse d'écoulement à l’intérieur de I’échantillon 
est provoquée par un gradient d’humidité. Par exemple, a l'instant 
t = 0, comme il n'existe de gradient d’humidité qu'à l'extrémité 
inférieure, au contact avec la membrane semi-perméable, la vitesse 
d'écoulement est nulle en chaque section, sauf pour s = 0. 

Dans cette méthode, on admet donc qu'il ne peut exister, à l’intérieur 
de l'échantillon, de gradient de potentiel s’il n'existe pas simultanément 
de gradient d’humidité. 

Dans le domaine d’expérimentation que nous avons exploré, c’est-à- 
dire aux faibles tensions, cette hypothèse ne semble pas en accord avec 
les phénomènes constatés car elle aboutit à des résultats aberrants. 
C’est ainsi que l'application systématique de la formule (11), à des échan- 
tillons de hauteurs différentes (1 ou 3 cm) aboutit à des valeurs du 
coefficient de conductivité capillaire fonction non seulement de la 
teneur en humidité mais de la hauteur de l'échantillon (cf. tableau IV 
et fig. 8). 

Il nous a paru nécessaire de revoir les hypothèses de calcul afin de 
déduire une autre interprétation des faits. 


Nouvelles hypothèses. 


Reprenons l'application de la loi de Darcy à l’ensemble de la phase 
liquide. Les limites spatiales de cette phase sont constituées par : 

— les ménisques, au contact de Pair 

— les particules solides, formant lignes de courant 

— la plaque poreuse, saturée. 


La fonction potentielle ou l’une de ses dérivées est également connue 
en ces limites. 


— elle possède une dérivée nulle, normalement à la paroi des parti- 
cules solides. 
— elle est égale à — w, à la partie inférieure de la plaque saturée. 


1 1 
— elle est égale à -ocosb (Œ + —) au niveau de chaque ménisque 
1 2 


de rayons de courbures principaux r, et ra, b étant l'angle de 
contact et o la tension superficielle du liquide. 
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Nous posons ` y = œ cosb d + by 
Ta Te 

y étant égal à y; à l'instant t = O et à y, lorsque tout écoulement 
a cessé. 

Si nous appliquons la théorie des écoulements à potentiel et la loi de 
Darcy définissant le coefficient de conductivité hydraulique à l’ensemble 
de la phase liquide, nous écrivons que la vitesse de drainage, au niveau 
de chaque ménisque est égale à : 


(32) 


dans laquelle : 


0 : teneur en humidité par unité de volume, en un point d’ab- 
scisse s 

: temps 

K : coefficient de conductivité hydraulique 

y : potentiel de capillarité correspondant aux rayons de courbures 
principaux du ménisque au point d’abscisse s (La fonction 
potentielle conditionnant l’écoulement est égale, en ce point, 
à — y) 

: direction de la ligne de courant issue du point d’altitude z 

U : surface normale aux lignes de courant, par unité de volume, 

au point d’altitude z 


+ 


EN 


En raison de l’équation de continuité, on a, en outre: 
dð = Uds (33) 


L'ensemble de deux relations (32) et (33) exprime l'application 
de la loi de Darcy à la phase liquide. 

D'autre part, en raison de la courbe d'équilibre (w, 0), nous con- 
naissons, pour chaque teneur en humidité, la valeur du potentiel de 
capillarité au niveau des ménisques. 

Nous admettons que cette relation reste valable lors du phénomène 
d'écoulement. 

Pour résoudre les équations d'écoulement, il faut pouvoir déterminer 


le gradient de potentiel à l’intérieur de la phase liquide. 
: Sg i z contovr nant 
En raison de la complication géométrique des lignes de courant es 


eemstitmest les particules solides, nous ne pouvons résoudre, algébri- 
quement, le problème de la répartition de la fonction potentielle à 
l'intérieur de l'écoulement. 
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Nous admettrons l'approximation suivante ` le gradient de potentiel 
dy i > e $ 
_ = est inversement proportionnel à z et directement proportionnel 
à la différence de la fonction potentielle existant aux deux extrémités 
de la ligne de courant, à l’intérieur de l’échantillon. 


Nous écrivons donc 


dy k ji 
V _ PTY on 


formule dans laquelle, outre les notations déjà définies, 


y  : valeur de la fonction potentielle, changée de signe, au point 
z = Ü d’intersection de la ligne de courant issue du point s et 
de la face supérieure de la plaque poreuse. 

z  : altitude comptée à partir de la face supérieure de la plaque 
poreuse. 


Nous admettons également que, dans l'intervalle de variation (y, 
vd on peut constamment écrire 
0-0 dé 
L= (35) 
Y—%Yo dy 
Yo et D correspondant à un point quelconque de la courbe (w, 0). 
Dès lors la combinaison des équations (32) (34) et (35) permet 
d'écrire : 


dé e o dy 
EF KU = Er (36) 
ou encore, en reprenant la notation de la formule (4), 
do 0-0 
= Fi z=0 
Er DU ` (36a) 


Les coefficients D et U varient en fonction de 9. Nous admettrons, 
comme précédemment, qu'on peut résoudre l’équation différentielle 
(36a) en utilisant une valeur moyenne DU, dans l'intervalle de varia- 
tion (0, 0,). 


Perméabilité infinie de la membrane. 


Si on utilise une membrane suffisamment perméable, on peut consi- 
dérer que son coefficient de Darcy K,, est infini. Dès lors le potentiel 
correspondant à z = 0 est constamment égal à — yp, et la teneur en 
humidité 0, _, devient immédiatement égale à 0, 


o 
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L'intégration de la relation (36a) fournit : 


z 0 0, 


t=— DU log, De? (37) 


la constante d'intégration ayant été calculée pour satisfaire aux condi- 
tions initiales et aux limites suivantes : 


t= 0 0<z<L 0 = 0, 
t = © Heset, 0 = 6, 
220 t>0 0 = 0, 


La quantité d'eau Q, qui traverse la membrane entre les instants 
t = O et t = t, est fournie par 


Q 
A 


PE L 
Lo, 0dz = L (0, —0)) (0— 0,) dz (38) 


CAL „ Ü 


dans laquelle la relation 0 (z, t) est fournie par la relation (37) ainsi 
écrite : 


9 0 Ra DUI 
(ES? =e * (37a) 
En combinant les relations (37a) et (38) et en posant 
J 
<= DUI (39) 
Z 
on aboutit à la relation : 
Qı DUt cs > 
1— > = e~ dy (40 
es “a WÉI 
Utilisant une notation classique (6, p. 867) nous posons : 
ES fv) = wo ez 
o fen = jus eg du (41) 
les fonctions Eo (v) et E, (v) étant connues par des tables (6). 
Dès lors, et si nous posons en outre, 
DU 
w =T (42) 
la relation (40) s'écrit : 
Qı p 
1— 3< = w [E(w) —E(w)] (40a) 
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Nous avons consigné, au tableau V, quelques valeurs de la fonction 
F(w) = w[Eo(w]) — Ex(w)] 
de façon à pouvoir en déduire une représentation graphique (figure 7). 


Calcul de DU. 
On construit à partir du tableau V un papier spécial (cf. fig. 8) 
dont l'échelle des ordonnées redresse la fonction (1 — 3) lorsqu’elle 


est exprimée en fonction de w. En raison de la relation (42), dans ce 
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diagramme, la fonction (1 — œ se traduira également par une droite 
f 


si on l'exprime en fonction du temps. Cette droite passe par les 
points 
Qi = 0, w = 0, t = 0 
L 


Q = 0,4Q„ w = 0,3854, t = 0,385 né 


Tantar V 


E (w) w [Eo(w)-E,(w) | 


0 — | — 1,0000 
0,01 99,00 4,038 0,9496 
0,02 49,01 3,355 0,9131 
0,03 32,35 2,959 0,8817 
0,04 24,02 2,681 0,8536 
0,05 19,02 2468 0,8276 
0,06 15,70 2,295 0,8043 
0,07 13,32 2,151 0,7818 
0,08 11,54 2,027 0,7610 
0,09 10,15 1,919 0,7408 
0,1 9,048 1,823 0,7225 
0,2 4.094 1,223 0,5742 
0,3 2,469 0,9057 0,4690 
0,4 1,676 0,7024 0,3894 
0,5 1,213 0,5598 0,3266 
0,6 0,9147 0,4544 0,2762 
0,7 0,7094 0,3738 0,2349 
0,8 0,5617 0.3106 0,2009 
09 04517 0.2602 0,1724 
1,0 0,3679 0,2194 0,1485 
1.1 0,3026 0,1860 0,1283 
1,2 0,2510 0,1584 O,1111 
1.3 0,2096 0,1355 0,0963 
1,4 0,1761 0,1162 0,0839 
1,5 0,1488 0,1000 0,0732 
1,6 0,1262 0,08631 0,0638 
17 0,1075 0,07465 0,0559 
1,8 0,09183 0,06471 0,0488 
1,9 0,07872 0.05620 0,0428 
2,0 0,06767 0,04890 0,0376 
3,0 0,01660 | 0,01305 0,0108 
4,0 0,004579 0,003779 | 0,0032 
50 | 0,001348 0,001148 0,0010 


10 


08- 
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Le report des points expérimentaux fournira un ensemble de points 
d'allure rectiligne. 

La droite de régression les représentant coupe la parallèle à l'axe des 
abscisses. 


Q: = 0,4 Q, 


en un point T qui permet de calculer DU par la relation : 


DU = 0,385 (43) 


10- 
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Nous avons reproduit (fig. 8) les courbes expérimentales déjà utilisées 
(fig. 3). : 

Nous avons également calculé le produit DU pour tous nos résultats 
expérimentaux (tableau IV). 


Comparaison des hypothèses émises. 


Rappelons la différence fondamentale. 

Dans la loi (10) qui découle de l'intégration de l'équation (4), on 
admet que le potentiel d’une particule d’eau, à l’intérieur de la phase 
liquide, est lié univoquement à la teneur en humidité moyenne du sol, 
au point considéré, ` 

Dans l'énoncé des nouvelles hypothèses, nous avons exprimé que la 
fonction potentielle ou ses dérivées ne sont connues, avec certitude, 
qu'aux limites de la phase liquide. 

Nous admettons ainsi, qu'à l’intérieur de cette phase, le potentiel d’une 
particule dépend des conditions aux limites et de la situation de cette 
particule à l’intérieur de la phase. 

D'autre part nous admettons que l'hypothèse du mouvement unidimen- 
sionnel ne permet pas une interprétation correcte des phénomènes. 

Moyennant cette restriction, on peut imaginer que les gradients 
locaux de potentiel ne correspondent pas, nécessairement, aux gradients 
locaux de la teneur en humidité. La relation (y, 0), traduisant la 
liaison entre la teneur en humidité et le potentiel (fig. 5) devient une 
loi d'équilibre hydrostatique. Elle n’est applicable, en équilibre hydro- 
dynamique, qu'aux seuls points des interfaces eau-air, c’est-à-dire au 
niveau des ménisques. 

Les conséquences de cette différence de conception sont les lois 
respectives exprimées par les relations (7) et (37a). 

La comparaison directe entre ces lois n’a pas été possible, dans le 
cadre de nos observations car nous n’avons pas mesuré, d’une façon 
continue, la variation de la teneur en humidité 9 au cours du temps. 

Force nous est de comparer les relations (10) et (40). 

Ces relations sont reportées graphiquement fig. 7. 

Le simple examen comparatif de cette figure et des diagrammes 
expérimentaux (fig. 3 ct fig. 8) montre qu'il est utopique d'espérer 
trancher entre les hypothèses conduisant respectivement aux formules 
(10) et (40) à partir de la forme des diagrammes expérimentaux seule- 
ment. 

Par contre, si nous exprimons les valeurs du coefficient DU calculé 
par la relation (43) et si nous exprimons ce coefficient en fonction de 0 
et de L (fig. 9), nous obtenons une dispersion des résultats beaucoup 
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moindre que dans le cas de D calculé par la relation (11b). Nous con- 
cluons en faveur des nouvelles hypothèses qui permettent une interpréta- 
tion des résultats expérimentaux plus cohérente et ne se heurte plus 
aux anomalies constatées précédemment. 


i 


IO o en cèl jour 


Évazuariox DE D ET DE K 


Dans l'écoulement envisagé, nousiestimons donc inapplicable l'équa- 
tion de l'écoulement unidimensionnel mise sous l'une des formes : 


50 
G == D 33 (8) 
60 
q = DS (8a) 
oy 
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Ceci ne signifie pas que la loi de DARCY soit infirmée. Au contraire en 
considérant le phénomène spatial il nous a été possible de définir un 
paramètre D ayant conservé la signification définie en fonction de la 
loi de Darcy (équation 4). 

Il apparaît intéressant de tenter une évaluation de D et de K, a 
partir de nos résultats. 

Ces paramètres ainsi obtenus pourraient être utilisés pour l'étude 
d’autres écoulements ou l'hypothèse du mouvement unidimensionnel 
peut être appliquée. 


ÉVALUATION DE U 


On ne peut en déduire une évaluation directe sans imaginer un 
modèle géométrique, toujours discutable ; la valeur de U ne sera donc 
connue qu'avec une certaine approximation. 

Nous utilisoñs le raisonnement, ci-dessous. 

Lors des mesures aux faibles tensions (y, < 40 cm) on constate que 
l'écoulement est négligeable. Le sol étudié reste, pratiquement, saturé. 

Dans la zone des écoulements importants (y, compris entre 50 et 
100 cm), il apparaît raisonnable d’admettre que les canalicules d’air 
qui s’infiltrent dans la phase liquide, peuvent être assimilés à un réseau 
de cavités de formes variables, reliées entre elles par des galeries cylin- 
driques. L'ensemble de cette phase gazeuse est tapissée par la phase 
liquide. Celle-ci, à son tour, enveloppe la phase solide. 

Les volumes respectifs de ces trois phases sont connus et sont reliés par 
la relation : 

0,=1-- 0- v (44) 
dans laquelle 

0, : volume d'air par unité de volume apparent de sol. 

H : teneur en humidité, exprimée en volume, par unité de volume 

apparent de sol. 

v : volume de la partie solide, rapportée à l'unité de volume 

apparent. 


On peut évaluer U, par excès, en admettant que la phase gazeuse est 
essentiellement formée d’un réseau de cavités identiques, deformes sphé- 
riques, le nombre de cavités étant égal au nombre de particules solides. 

Dans ce cas la surface spécifique U sera fournie par la relation : 


U = $ 023 vin 300 028 vin (45) 


dans laquelle d est le diamètre moyen des particules solides, exprimé 
en cm. 
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Cette relation n’est valable que si le rayon des cavités d'air est supé- 
rieur au rayon des ménisques correspondant a y. 

Admettons l’approximation suivante du rayon moyen r d’un mé- 
nisque 


r = = (46) 


On en déduit les valeurs numériques du tableau 6, les valeurs de 0 
utilisées étant celles de la courbe w de la figure 5 (Dessèchement). 


Tasreau VI 

w cm 30 | 3o | 70 90 110 

rem 0,00495 0,00297 0,00212 0,00165 0,00135 
9 cm*/cm* | 0.376 0,285 | 0,048 0,007 
Oa cm* cm* 0,004 0,109 | 0,319 0,379 

v cm?/cm? 0,620 0,606 0,633 0,614 
Went 5,12 579 | 120 | 134 

R cm | 0,00019 | 0,00 0,00738 | 0,00796 0,00852 


Ce tableau contient également les valeurs de U et du rayon R des 
cavités gazeuses correspondant à la formule (45). 

On voit que la restriction à faire intervenir lorsque r estsupérieur à R 
n’est pas à considérer dans le cadre des mesures rapportées (y > 50 cm). 

Il ne nous a pas été possible d'évaluer U par défaut. Nous ferons 
simplement remarquer que si le nombre de cavités gazeuses est dix fois 
moindre qu'il n’a été imaginé, la valeur de U doit être réduite dans la 
proportion de 1 à 2,15 = Vj. 

En conclusion, nous admettons les valeurs suivantes de U par inter- 
polation graphique des valeurs du tableau VI: 


yom 55 65 75 85 95 
Ucm 75 94 109 119 127 


On en déduit les valeurs du coefficient D renseignées à la dernière 
colonne du tableau IV. 


Conclusions pratiques quant à la méthode de mesure de D et de K 


Le tableau 7 et la figure 10 illustrent les divergences importantes 
auxquelles aboutissent l’utilisation de la méthode de GARDNER ou Pinter- 
prétation que nous proposons. (Le tableau 7 et la figure 10 fournissent 
les valeurs moyennes calculées à partir des résultats consignés au tableau 4). 
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"Tanta VII 


= Lem 


D U en Dencm| Ken Denempi Ken |D U en 
cm jour? cm jour cm cm 
jour- (méthode) ` Jour? |(méthode| jour! | Jour) 
| propo- Gardner) 


sée) 


25,7 f 044 | 0441.10 | 119,3 
32,3 0,31 | 0,30.10- | 125,8 
90,5 0,67 | 0,18.10-2 | 232,3 

164.9 | 1,29 | 0,39.10-2 | 3646 

205,7 1,84 | 0.42.10-2 | 428,6 


es TR UX D UL 
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